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ABSTRACT 

Binary stars in eccentric orbits are frequently reported to present 
increasing levels of activity around periastron passage. In this 
papeiu we present results of a calculation from first principles of 
the velocity field on the surface of a star that is perturbed by a 
binary companion. This allows us to follow the orbital phase- 
dependence of the amount of kinetic energy that may be dissi- 
pated through the viscous shear, E, driven by tidal interactions. 
For stars with relatively small stellar radii compared with the 
orbital separation at periastron (R\/r per < 0.14), a clear max- 
imum occurs before periastron passage for sub-synchronous ro- 
tation and after periastron for super-synchronous rotation. For 

larger stellar radii however, E oscillates over the orbital cycle 
and periastron passage does not cause a particularly greater en- 
hancement in energy dissipation rates than some of the other or- 
bital phases. Finally, we perform exploratory calculations for a 



1 Based on the poster presented at Hot Massive Stars... A lifetime of Influence, Work- 
shop held at Lowell Observatory, 2008 October 12-15, celebrating Peter Conti's birthday 
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WR/LBV binary system that in 1993-1994 underwent an LBV- 

like eruption, HD 5980. Our E calculations reproduce the os- 
cillations that appear around periastron passage in HD5980's 
recent visual light curve. We suggest that the energy dissipa- 
tion from tidal flows in asynchronously rotating binary stars may 
provide a mechanism by which evolving stars may be driven into 
an active state. Given the nature of the tidal perturbations, the 
resulting mass-loss distribution is expected to be non-uniform 
over the stellar surface and highly time-dependent. 



1 Introduction 

A number of binary systems show evidence of enhanced activity 
around periastron passage. Among these, T] Car is the most ex- 
treme and best documented example of periodic brightening at 
X-ray, visual and IR wavebands associated with periastron pas- 
sages. Recently, van Genderen & SterkerJ^ suggested that these 
periastron events may have the same physical cause as the milder 
"periastron effects" exhibited by many less renowned eccentric 
binaries in which small enhancement (Am„ ~0.01-0.03 ma5 ) in 
the visual brightness of the system around periastron passage 
are observed. They suggested that the fundamental cause of 
the effects may reside in the enhanced tidal force present during 
periastron passage. 

Our study of tidal interactions and their potential role in 
sparking stellar activity was first inspired by the peculiar be- 
havior of the WR/LBV system HD 5980@ In particular, the 
orbital phase-dependent behavior of the line-profile variations 
led us to raise the question of whether the periastron passage 



2 2007, IBVS, 5782 

3 Koenigsberger, G. 2004, RMAA,40, 107; Georgiev, Hillier & Koenigsberger, this meet- 
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Figure 1: Artistic rep- 
resentation of the one- 
layer model for the tidal 
interaction calculations. 
The inner body of the 
star is assumed to ro- 
tate uniformly, while the 
perturbations acting on 
the surface layer pro- 
duce a velocity field lead- 
ing to local angular ve- 
locities that differ from 
the underlying uniform 
rotation. 



could affect the intrinsic stellar wind properties, possibly en- 
hancing the mass-loss rate around this phaseE- Our objective 
is to explore whether the mechanism involved in causing the in- 
stabilities is the dissipation of kinetic energy of the viscous flows 
that are driven by the tidal interactions. 



2 Method for E calculation 

A system is in synchronous rotation when the orbital angular 
velocity Q equals the angular velocity of rotation uq. In eccen- 
tric orbits, the degree of synchronicity varies with orbital phase. 
We use periastron passage as the reference point for defining 
the synchronicity parameter (3 per = tuo/{} per = 0.02P v ™*^~^j 1/2 , 
where e is the orbital eccentricity, the rotation velocity v ro t is 

4 Koenigsberger et al. 2002, ASPCS 260, 507 
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given in km/s, the orbital period P is given in days, and the stel- 
lar equilibrium radius Ri is given in solar units. When at any 
orbital phase (3 ^1, non-radial oscillations are excited, driven 
by the tidal interactions. We refer to the azimuthal component 
of the forced oscillations as "tidal flows" . 

Our method consists of computing the motion of a Lagrangian 
grid of surface elements distributed along a series of parallels 
(i.e., rings with different polar angles) covering the surface of 
the star with mass Mi as it is perturbed by its companion of 
mass M2. The main stellar body below the perturbed layer is 
assumed to have uniform rotation. The equations of motion 
that are solved for the set of surface elements include the gravi- 
tational fields of Mi and M2, the Coriolis force, and gas pressure. 
The motions of all surface elements are coupled through the vis- 
cous stresses included in the equations of motion. The kinetic 
energy of the tides may be dissipated through viscous shear, 
thus leading to an energy dissipation rate, E. The rate of en- 
ergy dissipation per unit volume is given by the matrix product 
E = — P ?? : Vv, where P ?? is the viscous part of the stress tensor 
and v is the velocity of a volume element with respect to the 
center of the starj§. 

The benefits of our model are: 1) we make no a priori as- 
sumption regarding the mathematical formulation of the tidal 
flow structure since we derive the velocity field v from first prin- 
ciples; 2) the method is not limited to slow stellar rotation rates 
nor to small orbital eccentricities; and 3) it is computationally 
inexpensive. 



5 For details see Moreno & Koenigsberger 1999, RMAA, 35, 157; Moreno et al. 2005, 
A&A, 437, 641; and Toledano et al. 2007, A&A, 458, 413 
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Figure 2: Energy dissi- 
pation rates plotted 
function of orbital phase 
for a Mi =70 M Q pri- 
mary star of different 
stellar radii. In all 
cases, the star is in syn- 
chronous rotation at pe- 
riastron (<f) =1). E is 
normalized to its maxi- 
mum value over the or- 
bital cycle. P=19.3d and 
e=0.3. 
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Figure 3: E plotted as a 
function of orbital phase 
for a 70 M , Ri =10 
R binary system with 
stellar rotation velocities 
19-438 km/s. Maxi- 
mum E values always oc- 
cur within ±0.1 in phase 
from periastron. 



3 Dependence on stellar radius and rotation 
velocity 



For simplicity, we first consider a primary star that is syn- 
chronously rotating at periastron; i.e., (3 per =1. Figure 2 il- 
lustrates the behavior of energy dissipation rates as a function 
of orbital phase for a 70+54 M binary system with P =19.3 d, 
e =0.3 and values of Ri =10-22 R . For small radii, the equilib- 
rium tide component dominates. Thus, the trend for increasing 
E as periastron is approached reverses just before periastron, 
and goes to minimum at this phase since this is when (3 per =1. 
For stars with larger radii, however, the dynamical tide domi- 
nates, and E presents an oscillatory behavior over the orbital 
cycle. 

Figure 3 illustrates the result of holding constant Ri =1OR 
and varying v rot from sub-synchronous to highly super-synchronous 
values. The maximum E occurs close to periastron passage, but 
not necessarily centered at this orbital phase. We refer to the 
maximum value attained by E over the orbital cycle as Eyiax- 

Figures 4-6 summarize the behavior of Euax as a function 
of stellar radius and of rotation velocity for the set of mod- 
els listed in the Table, and Figure 7 plots the phase when this 
maximum is attained. Euax occurs before periastron passage in 
sub-synchronous starsr and after periastron when stars rotate 
super-synchronously. 
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Figure 4: Euax plotted 
aS cL function of stellar 
radius for the j3 per =1 
calculations listed in 
the Table. p is the 
mean density of the 
layer where the energy 
is dissipated. Squares 
correspond to cases 
with P=19.3, e=0.3 and 
triangles to P =100 d, 
e =0.767. 
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Figure 5: Same as pre- 
vious figure, but here 
the abscissa is given in 
units of stellar radius 
normalized to the or- 
bital separation at peri- 
astron. The mass-ratio, 
q =M 2 /Mi is listed. Col- 
ors are as in Fig. 4. 
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Figure 6: Emox for the 
binary system models 
with P=19.3d and e=0.3 
plotted function of 
equatorial rotation ve- 
locities. Note that ever- 
increasing rotation ve- 
locity does not produce 
an ever- increasing E Max . 
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Figure 7: Orbital phase 
when Euax occurs for 
the P per =1 cases for 
Mi = 10 M (red); 
40 M (blue); 70 M 
(green) . The dotted 
line shows that for stars 
with Ri/i peT < 0.14, 
EMax occurs before peri- 
astron passage while for 
stars with larger rela- 
tive radii, it occurs after 
periastron. 
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Table 1: Summary of binary system models 
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4 Periastron effect in HD 5980 ? 



Figure 8 is a comparison of the recent visual light curve of 
HD 59800 (kindly provided by P. Massey) with the E values 
obtained from our model calculation for a 70+54 M binary 
with (3 per =1.33 and Ri =21 R . The model qualitatively re- 
produces the behavior after periastron, suggesting that the tidal 
flows mechanism may be playing an important role in this sys- 
tem. Particularly interesting also is the apparent increase in 
the overall energy flux emitted around periastron ((f) =0.15) in 
the wavelength range AA1000-10000 A. Illustrated in Figure 9 
is the ratio of the spectral energy distribution from HST/STIS 
observations of 1999 and FUSE observations of 2002. This indi- 
cates that either there is a significant redistribution in the SED 
around periastron and/or there is an added source of energy. 
Hence, HD 5980 appears to belong in the list of binaries that 
present a periastron effect. Furthermore, we speculate that tidal 
flows may also provide a mechanism to explain the occasional 

6 Foellmi et al. 2008, RMAA, 44,3 
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Figure 8: V-magnitude 
light curve of HD 5980 
(from P. Massey) com- 
pared with the tidal en- 
ergy dissipation curve 
assuming =21 R 

and Po =1.33. Times 
of eclipses and periastron 
passage are indicated, as 
well as phases when the 
1999 HST/STIS spectra 
were obtained. 
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luminous blue variable-like eruptive behavior. 

We suggest the following working scenario: If the erupt- 
ing star in HD 5980 is in the process of leaving the Main Se- 
quence^, its systematically increasing radius should lead to val- 
ues of E that grow exponentially, as shown in Figs. 4 and 5. 
At some point, a critical radius could be reached at which the 
sum of E and the intrinsic luminosity of the star becomes super- 
Eddington, thus producing an episode of strongly enhanced mass- 
loss, which we would observe as an eruptive event. The fre- 
quency of this type of episodes might be expected to depend on 
the expansion rate of the star, the depth of the layer in which 
most of the energy dissipation takes place and the amount of 
mass that is ejected each time the critical radius is attained. 



7 as implied by its He/H abundance, Koenigsberger et al. 1998, ApJ, 499, 889. 
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Figure 9: Ratio of 
HST/STIS spectra ob- 
served at orbital phases 
0.15 and 0.83 in 1999. 
Error bars are rms of 
mean for wavelength 
range indicated by error 
bars for abscissa. 



5 Discussion and future directions 

The periastron events of rj Carinae, WR 140 and WR 125 are 
generally modeled in the context of interacting winds theory. § 
These models have been mostly successful, although some dis- 
crepancies remainj^ It is interesting to note that one of the dis- 
crepancies involves the time of maximum X-ray emission, which 
does not always appear at periastron passage as predicted by the 
wind-wind collision models. An explanation for this discrepancy 
may reside in that the mass-loss rate and wind velocity structure 
of the stars does not remain constant over the orbital cycle. As 
shown above, the tidal flows model predicts a significant pertur- 
bation of the stellar surface which, in turn, could significantly 

8 see Parkin & Pittard, 2008, MNRAS, 388,104, and references therein. 
9 see, for example, Rauw, G., 2008, RMAACS, 33, 59. 
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affect the stellar wind properties of the stars. 

A further consequence of the tidal flows model is that the 
strongest perturbation should occur on the larger of the two 
stars in the binary system. In the case of the WR+O binary 
systems, the O-star companion generally has the larger stellar 
radius. Hence, in a system such as WR 140 it is conceivable 
that periastron passage could lead to a significant increase in 
the O-star 's wind momentum. Assuming this to be the case, 
the possibility then exists that the O-star's wind could dominate 
the momentum ratio around periastron, thus even inverting the 
orientation of the shock cone geometry. 

Finally, it is interesting to note that the action of tidal flows is 
not distributed symmetrically over the stellar surface and thus, 
if the energy dissipated by these motions contributes towards 
mass-loss, the additional outflow is expected to depart from 
spherical symmetry as well. Furthermore, the critical radius 
could be reached by the star several times during its late stages 
of evolution, thereby producing multiple episodes of enhanced 
mass-loss prior to its demise in the supernova event. An inter- 
esting problem would be to analyze the morphology of circum- 
stellar structures that may be produced from such asymmetric 
mass-ejections as well as the manner in which this distribution 
of matter then may affect the expansion of the supernova ejecta. 

In conclusion, the tidal flows produced in asynchronously ro- 
tating binary systems provide a mechanism that may explain a 
wide variety of observational phenomena ranging from periodic 
photospheric line-profile variability to episodic and asymmetric 
mass-shedding events. The actual magnitude of the energy dis- 
sipation rates depends strongly on the viscosity parameter that 
is assumed for the stellar material and on the depth of the layer 
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where the energy dissipation rates are most significant. A proper 
assessment of these physical parameters will help determine the 
actual importance of the role being played by asynchronous ro- 
tation and tidal flows in producing peculiar phenomena observed 
in binary stars. 
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